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Advanced  (non-simple-dipole)  EMI  models’  inversion  and  classification  performance  is 
presented  for  the  ESTCP  Live-site  UXO  Discrimination  Study  at  former  Camp  Butner,  NC. 
The  advanced  models  combine:  (1)  the  joint  diagonalization  (JD)  algorithm  for  estimating 
number  of  potential  anomalies  from  the  measured  data  without  inversion,  (2)  the  ortho- 
normalized  volume  magnetic  source  (ONVMS)  for  representing  targets’  EMI  responses  and 
extracting  targets’  intrinsic  parameters  feature  vector,  and  (3)  the  Gaussian  Mixture  algorithm 
and  probability  neural  network,  that  utilizes  the  extracted  discrimination  features  for  classifying 
buried  objects  as  targets  of  interest  or  not.  Namely,  the  studies  were  conducted  for  the  next 
generation  sensor  data:  Time-domain  Electromagnetic  Multi-sensor  Towed  Array  Detection 
System  (TEMTADS)  and  Metal  Mapper  (MM)  sensors’  cued  data  sets  collected  at  the  Camp 
Bunter,  live  UXO  site.  These  sensors  provide  the  measured  multi-static  response  (MSR)  data 
matrix.  Eigenvalues  versus  time,  which  are  determined  using  the  JD  from  the  MSR  data  matrix, 
provide  information  about  the  number  of  targets  contributing  to  the  signal  and  their  initial 
classification  features.  Once  the  number  of  targets  is  known,  then  data  are  inverted  and  intrinsic 
parameters,  such  as  the  total  ONVMS  that  is  a  function  of  target’s  geometry  and  material 
composition,  are  determined  for  each  potential  target.  These  intrinsic  parameters  are  grouped 
using  the  unsupervised  Gaussian  mixture  approach.  For  each  group  an  anomaly  is  identified  and 
ground  truth  is  requested.  Once  the  requested  ground  truth  data  are  obtained,  then  each  of  the 
groups  is  classified.  In  this  presentation,  the  advanced  EMI  methods’  data  inversion,  processing 
and  discrimination  scheme  will  be  reviewed,  and  the  classification  results  scored  by  the  Institute 
for  Defense  Analyses  (IDA)  will  be  presented  for  both  the  TEMTADS  and  MM  sensors  Camp 
Butner,  NC  cued  data  sets. 
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>  Advanced  EMI  Models 

□  Normalized  Surface  Magnetic  Source  (NSMS)  Model 

□  Orthonormalized  Volume  Magnetic  Source  (ONVMS)  Model 

>  EMI  Data  Pre-processing  and  Inversion 

□  Joint  Diagonalization 

□  Direct  Search  technique  for  Multi  Targets 

>  Classification 

□  Clustering 

□  Library  Matching 

□  Results  (IDA  Score) 

>  Summary 
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UXO  classification 


The  entire  UXO  classification  process  can  be  divided  into  three  parts 


1.  Data  Collection 


2.  Data  Inversion 


3.  Decision 


Forward  Models:  NSMS 


Normalized  Surface  Magnetic  Source  Model 


Primary  field  induces  eddy 
currents  inside  metallic  objects 


>  Extended  dipole  model 


>  NSMS  model  accounts  for  target’s  heterogeneity. 

>  Total  NSMS  is  an  intrinsic  target  parameter. 
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NSMS  Applied  to: 


1.  APG  test  site  (214  anomalies) 

APG  Discrimination  results  were  excellent: 

>  All  UXO  items  were  correctly  identified  as  TOI 

>  All  TOI  items  were  correctly  identified  by 
type/caliber 

>  There  was  a  5  %  false  positive  rate 


APG  Targets  of  interest  (TOI) 


1 1 1 

25- mm  37  -  m  m  60  -  m  m 


81 -mm  105 -mm  105- mm  HR 


2.  SLO  Live  UXO  site 


SLO  Discrimination  results:  SL0  Tar9ets  of  interest  (TOI) 


>  One  false  negative  for  Metal  Mapper  (2492  anomalies). 

>  Seven  false  negatives  for  TEMTADS  (1464  anomalies). 


The  main  challenges  there  were: 

>  Multiple  overlapping  targets 

>  low  signal  to  noise  ratio 


60-mm  2.36  inch  81  mm  4.2  inch 
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Forward  Models:  ONMVS 

Ortho-Normalized  Volume  Magnetic  Source  (ONVMS)  Model 


The  ONMVS  model  divides 
the  computational  space  into 
cells. 


The  key  elements  of  the 
ONMVS  are: 


>  The  scattered  EMI  field  is  approximated  using  an  orthonormalized 
function  expansion: 

_  r ,  =T  =  \  0 ,  i*k 

H(r)  =  ^<//,(r)V>;,  Where  ){Yi  ■Yk)dv  =  \ 

F.,m  =  k 


i=\ 


ONVMS  continued 

>  These  orthogonal  functions  are  constructed  using  the  scattered 
magnetic  field’s  Green  function  via  Gram-Schmidt  ortho-normalization 
process: 

r)  =  <T,(r)-  X  where  for  i  <  k,  Alk  =  0, 

k  =  1 

>  The  modeled  Magnetic  field  is  fitted  to  measured  data,  and 

>  The  targets  features  are  extracted 


The  ONVMS: 

□  Avoids  an  ill-conditioned  matrix; 

□  Separates  overlapping  targets  easily; 

□  Provides  total/effective  polarizabilities; 

□  is  applicable  for  non-uniform  sub-volumes. 
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Joint  Diagonalization 


Time  Channels  1:121  Representing  TEMTADS  data  in  Space  -Time 


Hd  =VTD(tn)VT 


Hd  =VD(t2)VT 


t 

H,  -VD(t,)VT  | 

*  V  is  a  Eigen  | 
vectors 
matrix 


Tx:  1:25 


•  D(t)  is  a 
diagonal 
matrix 


^  Rx:  1:25 


The  goal  is  to: 

>  determine  the  eigenvalues  of  Hd  tensor  for  each  time  channel. 

>  find  an  eigenvector  V  that  will  be  shared  by  all  matrices. 


D(tk)  =  VTHd(tk)V,  k=l,  2, n 
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Eigen  Values  [Arb] 


JD  applied  to  CB-TEMTADS  data; 

eigenvalues  for  classification 
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The  eigenvalues  show  targets  features: 

Two  targets  (105  mm  HE  and  105  mm 
HEAT),  having  same  size,  but  different 
material  properties  have  different  time 
decaying  eigenvalues 


Eigen  Values  [Arb] 


JD  applied  to  CB-TEMTADS  data; 

eigenvalues  for  classification 


TOIs  have  slow  time 
decaying  eigenvalues 

37  mm  projectiles  with 
copper  bands  have 
distinguishable 
eigenvalues 
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Eigen  Values  [Arb] 


JD  applied  to  CB-TEMTADS  data; 


eigenvalues  for  classification 


TOIs  have  slow  time 
decaying  eigenvalues 

Clutter  items  have  fast 
time  decaying 
eigenvalues 


Eigen  Values  [Arb] 


JD  applied  to  CB-TEMTADS  data; 

eigenvalues  for  Multi  targets 
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Eigenvalues  [Arb] 


JD  applied  to  CB-TEMTADS  data: 

Resolving  small  signal  to  noise  ratio 
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*  eigenvalues  above  threshold 


10'1  10°  101 
Time  [sec] 


The  eigenvalues  are  small,  but  decay  slowly  in  time,  that  means: 
the  anomaly  is  buried  deep  and  it  is  a  potential  TOI. 
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CB-TEMTADS  Data  Classification  Approach 


>  JD  applied  to  all  2293  CB-TEMTADS  data; 

>  The  number  of  potential  targets  were  estimated  using  JD; 

>  The  first  Dig  list  was  created  based  on  Eigenvalues, 

>  All  data  sets  were  inverted  using  the  ONVMS  technique; 

>  The  effective  polarazabilities  were  determined 

>  70  Custom  training  data  sets  were  requested; 

>  Targets  were  ranked  via  Library  matching; 


Total  ONVMS  [Arb] 


Library  Matching  applied 


to  CB-TEMTADS  data 


Time  [msec] 


Time  [msec] 


Time  [msec] 
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Camp  Butner  TEMTADS  Classification  results 


Scored  Results  for  the  TEMTADS  Cued  Data  Sets: 


Butner  Dartmouth  AdvancedModels  None  TEMTADS  Custom  via  TOI 


Difficult  TOIs: 
#460  (37mm) 
#2504  (37mm) 
#2405  (37mm) 
#2340  (37mm) 
#2023  (37mm) 
#2017  (37mm) 
#1981  (37mm) 
#1945  (37mm) 


Butner  Dartmouth  AdvancedModels  None  TEMTADS  Custom  via  TOI 


Number  of  Non-TOI  Digs 


>A!I  data  were  inverted  and 
analyzed. 

>  No  False  Negatives:  all 
TOI  were  indentified 
correctly. 

>AI1 105  mm  and  37  mm 
were  identified  by 
caliber/type; 
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CB-Metal  Mapper  Data  Classification  Approach 


^  All  data  sets  were  inverted  as 

□  One 

□  Two 

□  Three 

targets  using  the  ONVMS, 

^  The  effective  polarazabilities  were  determined; 


^  Targets  were  clustered  using  the  principal  effective  polarazabilities  ; 


Pol  [Art] 


ONVMS  applied  to  CB-Metal  Mapper 

data:  Anomaly  #2504 


Single  source  inversion 
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Multi  targets  inversion 


18 


Pol  [Art] 
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ONVMS  applied  to  CB-Metal  Mapper 

data:  Anomaly  #2405 


Single  source  inversion 
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■  Case#2405oGSo  1  oGeo 


37  mm:  Library 


Inverted  parameters 
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Multi  targets  inversion 


■Case#2405oGSo1oGeo 


37  mm:  Library 


Inverted  parameters 
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Gaussian  mixture  model  for  MM-ONVMS 

clustering 


>  Uses  discrimination  features  from  ONVMS 

>  Builds  the  mixture  Gaussian  distribution  for  K  clusters; 

>  The  expectation-maximization  algorithm  used  to  estimate 
weight,  mean  and  variance  for  each  of  the  K  clusters; 


le 


105  mm -HEAT 


37mm  -2 


105  mm  -HE 


37mm  -1 


Mzz(t,)  and  Mzz 
(*i  yMzz(t30)  are 
used  as 
discrimination 
features. 


INERT" 
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Camp  Butner  Metal  Mapper  Classification  results 

Scored  Results  for  the  Metal  Mapper  Cued  Data  Sets: 


Butner  Dartmouth  Advanced  Models  None  MetalMapper  Custom  v2  TO  I 


Difficult  TOIs: 
#884  (37mm) 
#178  (37mm) 
#1728  (Fuze) 
#1002  (37mm) 
#1000 (Fuze) 
#856  (37mm) 
#773  (37mm) 
#720  (37mm) 


Butner  Dartmouth  Advanced  Models  None  MetalMapper  Custom  v2  TO  I 


Number  of  Non-TOI  Digs 


>121  Custom  training 
data  sets  requested 

>AII  data  were  inverted 
and  analyzed. 

>  No  False  Negatives. 

>AII  TOI-s  were  identified 
by  caliber/type; 
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Summary 


>  Advanced  EMI  models  applied  to  CB  Cued  Data  sets. 

>  The  Models  are  robust  and  noise  tolerant. 

>  They  are  applicable  for  single  and  multi  targets. 

>  Classifications  are  done  using  LM,  JD  and  Gaussian 
mixture  clustering. 

>  Excellent  classifications  were  demonstrated. 

>  No  False  Alarms. 

>  The  models  are  adapted  for  all  advanced  EMI  sensors. 

>  The  technology  will  be  tested  further  under  the  new 
ESTCP#  MR-201 101. 
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